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ABSTRACT 

Recent measurements of Si iv/C IV ratios in the high-redshift Ly a forest [Songaila & Cowie, AJ, 1 12, 335 
(1996a); Savaglio et al , A&A (in press) (1997)] have opened a new window on chemical enrichment and 
the first generations of stars. However, the derivation of accurate Si/C abundances tequires reliable 
ionization corrections, which are strongly dependent on the spectral shape of the metagalactic ionizing 
background and on the “local effects” of hot stars in nearby galaxies. Recent models have assumed 
power-law quasar ionizing backgrounds plus a decrement at 4 Ryd to account for He n attenuation in 
intervening clouds. However, we show that realistic ionizing backgrounds based on cosmological radiative 
transfer models produce more complex ionizing spectra between 1 -5 Ryd that are critical to interpreting 
ions of Si and C. We also make a preliminary investigation of the effects of He II ionization front 
nonoverlap. Because the attenuation and reemission by intervening clouds enhance Si iv relative to C IV. the 
observed high Si iv/C iv ratios do not require an unrealistic Si overproduction [Si/C>3 (Si/C)©]. It the 
ionizing spectrum is dominated by “local effects” from massive stars, even larger Si iv/C iv ratios are 
possible. However, unless stellar radiation dominates quasars by more than a factor of 10, we confirm the 
evidence for some Si overproduction by massive stars; values Si/C^2 (Si/C)© fit the measurements better 
than solar abundances. Ultimately, an adequate interpretation of the ratios of C iv, Si iv, and C il may 
require hot, collisionally ionized gas in a multiphase medium. © 1997 American Astronomical Society. 
[S0004-6256(97)00205-7] 


I. INTRODUCTION 

The detection of metals associated with the “high” col- 
umn density [N{H I) > 10 l45 cm -2 ] Lya forest of absorbers 
in quasar spectra (Cowie et al 1995; Tytler et al 1995) pro- 
vides new challenges to the understanding of the nature and 
origin of these absorbers. The contamination by metals must 
now be accounted for in models for cloud evolution. It is 
also possible to perform a detailed analysis of the thermody- 
namic and ionization state of these absorbers, which have 
previously been probed only by their hydrogen Lya absorp- 
tion and the integrated effect of their He n Lya absorption 
(Giroux et al 1995). If the clouds are photoionized, the de- 
tection of C II, C IV, and Si IV may constrain the shape of the 
metagalactic radiation spectrum and indicate the epoch ot 
dramatic changes in the ionizing radiation field at z > 3 (Son- 
gaila & Cowie 1996a; Savaglio et al 1997). 

The multiple C iv lines associated with a single H I ab- 
sorption complex (Songaila & Cowie 1996a) belies a simple 
picture of a homogeneous cloud. These absorbers, like Ly- 
man limit systems, may better be interpreted with multiphase 
models (Giroux et al 1994; Petitjean et al 1996; Haehnelt 
et al 1997; Hellsten et al 1997) that may include hot, colli- 
sionally ionized gas. We defer such interpretation to a later 
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study, focusing here on the tradition of analyzing metal-line 
ratios with single-phase photoionization models. In this way, 
we can isolate the effects of the incident spectrum on species 
with ionization potentials in the range 1-5 Ryd (e.g., C in. 
C IV, Si ill. Si iv). For example, Songaila & Cow ie (1996a) 
argue that the observed increase in the Si iv/C IV ratio at 
£>3.1 suggests that the Hem regions around ionizing 
sources may not have overlapped by that epoch, cutting off 
all radiation with h v > 54.4 eV. Here, we primarily consider 
models for the averaged metagalactic background in which 
He 11 ionization fronts have overlapped before We do 

consider a limiting case in which no photons above 4 Ryd 
are present in the background, as well as a case where we 
restore the much less heavily attenuated x-ray background. 

In the past decade, metal-line ratios in Lyman Limit 
(A^hi— 10 17 cm 2 ) absorbers have been used to probe the 
shape of the ionizing spectrum, and to estimate the intrinsic 
metallicity of the absorbers (cf. Bergeron & Stasinska 1986; 
Sargent et al 1988; Steidel et al 1988; Steidel & Sargent 
1989; Bergeron & Ikeuchi 1990; Miralda-Fscude Sc Oslriker 
1990; Donahue & Shull 1991 ; Madau 1992). With the advent 
of new spectrographs on large telescopes, the lower N H , ab- 
sorbers can be exploited in the same way. These absorbers 
hold a special interest, because the source of their metal en- 
richment is still unclear. While the metal lines in Lyman 
limit absorbers have been associated with the haloes of 
bright galaxies, the number of Lya forest absorbers greatly 
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exceeds the corresponding number of bright galaxies, and the 
Lya absorbers do not seem to show the same clustering 
properties as galaxies. 

As yet, metal lines have only been associated with “high- 
column’' absorbers [N( H l)> 10 14 5 cm 2 ], which numeri- 
cally are a small fraction of the Lya forest, it is still possible 
that these clouds may represent an overlap of the high- 
column end of the pristine Lya forest clouds and the low- 
column end of the metal enriched clouds associated with 
galactic haloes. As a result some, if not all, of the metal-line 
absorbers may still be associated with the galaxies that con- 
tain the stars responsible for their metals. For example, over 
a Hubble time at z = 3, galactic outflows at 300 km s~ 1 
will transport heavy elements a distance 

./ 1 +z\~ 3a 

d*(250 kpc^Wr 1 )— ) (1) 

for a Friedmann universe with H 0 =( 100 km s _ l )h and Cl Q 
= 1 . In practice, the period of heavy element injection may 
have lasted only 10 8 yr, and the metals could move dis- 
tances of only 30-50 kpc from their sources (bright galaxies, 
dwarf galaxies, globular clusters). 

As Madau & Shull (1996) have emphasized, metal enrich- 
ment of the Lya absorbers implies a substantial Lyman con- 
tinuum (LyC) emission accompanying the star formation at 

3.5. Thus, a large fraction of the metagalactic background 
at high redshift could be due to massive stars, a point also 
made by Giroux & Shapiro (1996) and Savaglio et al 
(1997). Furthermore, the absorbers with metal lines may be 
sufficiently close to the local sources that the incident radia- 
tion field is dominated by local hot-star radiation. In contrast, 
if the metal enrichment of the Lya forest arises from Popu- 
lation III stars at z> 10 (Couchman 1985; Ostriker & Gnedin 
1996), the stellar ionizing radiation will be greatly attenuated 
by intergalactic absorption, and the metagalactic background 
at z — 3 will be dominated by the harder spectrum of quasars. 
Because metal-line systems would then be less likely to be 
associated with nearby sources of stellar radiation, their 
absorption-line ratios would be more representative of the 
metagalactic ionizing background. 

Preliminary analyses (Songaila & Cowie 1996a; Savaglio 
et al 1997; Haehnelt 1997) of the Lya forest clouds have 
assumed photoionization by simple power-law ionizing spec- 
tra. The only way these models account for filtering of the 
ionizing source spectrum by intervening clouds is to include 
a decrement at r HeII (4 Ryd). More realistic models include 
several additional effects. First, when absorption due to in- 
tervening clouds is included, the background spectrum just 
above the H l and He ll ionization edges shows a flatter 
power law than that of the sources (Miralda-Escude & Os- 
triker 1990; Madau & Meiksin 1994; Giroux et al. 1995). 
Second, Haardt & Madau (1996) and Fardal et al (1997) 
show, using cosmological radiative transfer, that cloud emis- 
sion by recombination to the H 1 and He H ground states and 
He il Lya and two-photon radiation alters the extent of the 
decrement at r Hcll and changes the shape of the spectrum 
between I and 5 Ryd. Including a more realistic spectral 
shape will alter the relative populations of metal ions with 
edges in the 1-5 Ryd range even if the ionizing sources are 


all assumed to be quasars with intrinsic power-law ionizing 
spectra. 

In this paper, we consider the extent to which these alter- 
ations of the assumed shape of the incident radiation field 
affect the interpretation of the metallicity of these absorbers 
and the nature of the photoionizing sources. We use the data 
sets of Songaila & Cowie (1996a) and Savaglio et al (1997), 
which include Si IV and C IV lines associated with higher 
column density Lya forest clouds. Songaila & Cowie 
(1996a) provide an important subsample for which C n lines 
are also observed. 

2. THE METAGALACTIC RADIATION BACKGROUND 
2.1 Ionizing Sources 

Assumptions about the nature of the ionizing sources of 
the metagalactic background are uncertain, since the intrinsic 
EUV spectra of quasars are not well defined. For simplicity, 
we assume that all quasars have constant power-law spectra 
F t , oc „-<**. Our choice of spectral index, a s = 1.8, is consis- 
tent with HST observations of the EUV spectrum of quasars 
(Zheng et al. 1997) for radio-quiet quasars at z— 1 .5-2.0. As 
discussed by Fardal et al (1997), if the ionizing background 
(1-5 Ryd) is dominated by these quasars, the He n Lya ab- 
sorption optical depth r HeiI ^ 1 at z = 3. For clarity, we de- 
note the cloud-filtered, quasar-dominated mean intensity by 

J Q,f( v )* 

The characteristic spectrum of stellar sources in starburst 
galaxies is also not well known. Most models agree (e.g., 
Bruzual & Chariot 1993) that few photons are present with 
frequency above 54.4 eV, and that the spectrum falls off 
above h v 45 eV. Recent line-blanketed, non-LTE hot-star 
models (Gabler et al. 1992; Sellmaier et al 1996; Schaerer 
& de Koter 1996) vary considerably in the amount of radia- 
tion in the He I continuum (24.6 < hv < 54.4 eV) they pre- 
dict, compared to the LTE line-blanketed models of Kurucz 
(1992). An ongoing uncertainty, for our purposes, is the lack 
of low-metallicity stellar models, which are more appropriate 
for high-redshift starburst galaxies. For the present work, we 
use a starburst spectrum (Sutherland & Shull 1997) dis- 
cussed in the context of the ionizing spectrum of metal- 
forming galaxies by Madau & Shull (1996) (see their Fig. 3). 
In brief, it is a time integrated spectrum of a Gaussian star- 
burst of width 2 a = 4 Myr containing 5000 . of stars be- 
tween 8 and 85 . with IMF slope r=1.6, where dN 
( > . /t)ld. M sc . /// ] and N( > . M) is the number of stars 
with mass greater than . //,. 

Estimates of the mean metagalactic ionizing intensity at 
z > 2, based on the statistical diminution in the number of 
Lya forest lines close to quasars (the proximity effect), indi- 
cate a specific intensity at the hydrogen Lyman limit, v H , of 
y _ 2 i — 0.7 — f .5, where J__ 2 \ is in units of 

10‘ 21 ergs cm -2 s~ l Hz~ l sr 1 (Cooke et al 1997). If all 
Lya absorbers possess a metallicity of Z=0.01Z o by red- 
shift z ^ 3.5, the stellar ionizing photons associated with the 
production of the metals may be comparable to the ionizing 
photons produced by quasars (Madau & Shull 1996; Giroux 
& Shapiro 1996). If these ionizing photons are produced at 
higher redshift (e.g.. Population III stars at z > 5) they will be 
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highly attenuated by IGM absorption at z > 4. A further un- 
certainty is the fraction, (/ esc ), of ionizing photons that suc- 
cessfully escape from the galaxy in which the stars reside 
(Dove & Shull 1994). It does seem possible, though, that 
stellar sources could be an important contributor to the meta- 
galactic background, which could dominate in regions local 
to starburst galaxies. 

For example, the total Lyman continuum production rate 
of the Milky Way was recently inferred from COBE obser- 
vations of [Nil] 205 (Jim emission to be 3.5 X 10 53 
photons s" 1 (Bennett et al 1994), and most starburst galax- 
ies exceed the value S = ( 10 54 )S 54 photons s' 1 . Assuming 
that a substantial fraction, (/ esc ) ^ 0.1, of these photons es- 
cape these starburst galaxies, the mean photoionization rate 
of hydrogen due to the galactic source may exceed J 2 i 
= 1 if the absorber lies within a distance 

R a =( 88 kpc)S ] £jZ'J(f e J l/2 j -Jj-j • (2) 

Here, we have assumed the QSO spectra in the range 1-5 
Ryd to be power laws with spectral index a s ** 1.8. 

In this paper, we account for a stellar contribution to the 
mean intensity falling upon our absorbers in a global way. 
We assume that up to 2/3 of the ionizing emissivity from 
primary sources is stellar, as opposed to secondary sources 
due to recombination radiation from clouds. The remainder 
of the primary emissivity is then from AGN. We then com- 
pute the radiative transfer of primary and secondary sources 
through the IGM to derive the metagalactic ionizing back- 
ground, Jsqj(v)- ^e mean intensities Jqj (pure QSO 
sources) and J SQ f (stars + QSOs) therefore bracket our esti- 
mates of the contribution of stellar sources to the uniform 
metagalactic background, J MB (>). However, even if the 
metagalactic ionizing background is quasar dominated, the 
absorber may primarily see a local unfiltered source of 
stellar photons, We consider models which assume 

•/|V*/mb=5, 10, 20 at v = lor J qj und J[JJ mb ^ 

5, 10, 20 for J mb = J sqj '■ 


2.2 He II Ionization Fronts 


Within a picture of ionizing sources turning on in a neu- 
tral medium, Miralda-Escude & Rees (1994) made the im- 
portant suggestion that the He n I-front which propagates 
away from sources with a sufficiently soft ionizing spectrum 
can lag the hydrogen I-front. Madau & Meiksin (1994) 
showed that, if a v = 1 .9, the overlap of He n I-fronts could be 
delayed to as late as z = 3. As a result, both Jakobsen et al 
(1994) and Songaila & Cowie (1996a) suggested that, at 
z=3.1. Hell I-fronts around. the dominant ionizing sources 
may not have overlapped. Although a full treatment of over- 
lapping He II I-fronts is beyond the scope of this paper, we 
can make a simple estimate of the ratio of I-front velocities 
for H I and He II. From the flux-limiting equations governing 
I-front propagation (Shapiro & Giroux 1987; Donahue & 
Shull 1987), the velocity relative to the local Hubble flow is 
given by 



47 Trjrii ' 


(3) 


where / refers either to the H I ionizing front (1 Ryd con- 
tinuum) or to the Hell ionizing front (4 Ryd continuum). 
(The He 1 ionizing front is assumed to be coincident w ith the 
H I front, and we neglect recombinations and attenuation 
within the ionized zone.) Here, S, (photons s' *) represents 
the photon production rates of the source in the H I and He II 
continua, n t is the density of H I or He II, and r { is the front 
distance measured from the source. 

The propagation speeds of the H I and He II ionization 
fronts will be governed by the “flux-to-density” ratio, 
{SJnj *•), at 1 and 4 Ryd. For a QSO spectrum with power- 
law index a s , one finds that S Ucn /S Hl = 4 ~ The primor- 
dial helium abundance by mass has been estimated at values 
Y P = 0.231 ± 0.006 (Skillman & Kennicutt 1993) and 
Y P = 0.232 ±0.003 (Olive & Steigman 1995), while Copi 
et al. (1994) have suggested a 2<x concordance range of Y P 
— 0.221 — 0.243 from theoretical models of Big-Bang nu- 
cleosynthesis. If helium has a cosmological abundance, 
n He In h — 0.0785 by number (Y = 0.239 by mass), the flux-to- 
density ratios at a fixed radius r t are equal for 


In 4 


1.84. 


(4) 


At this critical spectral index, the H I and He ll ionization 
fronts will coincide. Interestingly, this critical index is very 
close to the mean index found by Zheng et al. (1997) for 
radio-quiet quasars. For AGN wdth harder spectra (a s 
< 1.84), the He II ionization front will precede the H I front 
by a small amount, equivalent to a few optical depths in the 
4 Ryd continuum. However, for those AGN w ith softer spec- 
tra (a v >1.84), the Hell front whll propagate at a lower 
speed and lie well within the H I front. For example, Zheng 
et ai (1997) found that radio-loud quasars had a mean index 
a s *** 2.2. In this case, the regions outside the He II fronts will 
see no 4 Ryd continuum radiation. 

We therefore consider a limiting case in which all photons 
above 4 Ryd are removed from the metagalactic background 
spectrum. We also consider an intermediate case, arising 
from the fact that quasars emit radiation at least into the 
x-ray range of energies, so that high-energy photons will not 
be strongly attenuated, even in regions where helium is en- 
tirely in the He II ionization stage. In the Lya clouds, we 
estimate the effect of the presence of these high-energy pho- 
tons by presuming that H II regions have overlapped by z 
= 4, and that the Lya forest clouds, whose distribution in 
N h j is known (see Sec. 2.2). have a maximal He Ii/H ! ratio. 
In interpreting He il Gunn-Peterson measurements, several 
groups (Madau & Meiksin 1994; Giroux et al 1995) find 
that the ratio, rj = N( He U)/N{ H I), must he 50—100 in order 
to explain the X < 304 A absorption as He ll line opacity in 
the H l Lya forest. This large ratio arises because He II is 
more difficult to photoionize than H i. 

The H neutral fraction / H , is always less than 
(n He /n H )r/ \ where rj = AT He ll)/A(H I) w'ithin an absorb- 
ing cloud. Thus, the maximal ratio rj nydS is bounded by limits 
on the baryon density from Big Bang nucleosynthesis 
(BBN), il h h 2 0.024 (Copi et al. 1994). At z = 3.4, if only 
absorbers with /V H t < 10 !> cm 2 are considered to contribute 
to 0 Hl , we find tfmax < 5° ()0 - This ,imit produces a local 
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Table: I . Significant multistage metals. 


Species 

£,h(Ryd) 

C iun.iv 

24.38, 47.89, 64.49 

N i.iuiuv 

14.53, 29.60, 47.45, 77.47 

o uuu 

13.62, 35.12, 54.94 

Ne i, ii. m 

21.57, 40.96, 63.46 

Si ii,iii,iv, v 

16.35, 33.49, 45.14, 166.77 

Al H.m 

18.83, 28.45 

S luiuv.v 

23.33, 34.83, 47.3, 72.68 

Fe iuiuv 

16.19, 30.65, 54.8 


continuum optical depth at 4 Ryd, drldz > 180 at z = 3.4, 
falling to 4 by v = 25 j^ h , the same continuum optical depth 
as that at 1 Ryd. As a result, for this model with nonover- 
lapping He II I-fronts, we only remove photons between 4 
and 25 Ryd, retaining metagalactic radiation background for 
v > 25 . This is a more conservative estimate than if all He 

is assumed to be in He II and distributed uniformly. 


Table 2. Absorber distribulion. 


Range in N n 

A, 

a 

10 5 - 10 3 

0.132 

1.41 

10 ^-0.02 

6.31 X 10 3 

1.85 

0.02- I0 5 

2.9X 10‘ 2 

1.46 


Lya (X304) line radiation. Accounting properly for the Lyc 
and Bac radiation is critical to the determination of the level 
of the mean intensity above 4 Ryd. While the line and 
2-photon radiation are not important for that, they strongly 
alter the spectrum between 1 and 3 Ryd and affect the popu- 
lations of ions with thresholds in that range (see Table 1). A 
complete discussion of the treatment of reemitted radiation 
used in computing the filtered spectra is given in Fardal et al 
(1997). 

3. Lya FOREST CLOUD MODELS 


2.3 Cloud Opacity/Emission 


For a random spatial distribution, the local continuum op- 
tical depth of the intervening absorbers is given by (Paresce 
etal 1980) 


dr{ v) 
dz 


-r 

Jo 


d l N c 

dN H ,dz 


[\-e A, Hi< r efK>')]^ Hi> 


(5) 


where, for clouds composed of H and He, 


o- eff ( v) = a H ,( v) + 


N 


He [ 


N 


H i 


0Hei( V) + 


N 


He ii 


N 


H 1 


O'He llM* 


( 6 ) 

and d 2 N c ldN H] dz is the distribution of absorbers in column 
density and redshift. This distribution has often been param- 
eterized by a function of the form 


d 1 N c 
dNft \dz 


=A(i+zVN H P. 


(7) 


for N/(H I) ^ A Hi (cm -2 ) A W (H I). Using new line lists 
from Keck data, Fardal etal. (1997) have completed a re- 
view of the statistics available on Lya absorbers, Lyman 
Limit Systems, and damped Lya systems. They suggest that 
the distribution of absorbers is best represented by a broken 
power law in N H , . If 


d 2 N c 
dN ]7 dz 


=A i (i+zVN;f‘, 


( 8 ) 


where A 17 = N H ,/( 10 17 cm -2 ), y — 2.585, and and A, rep- 
resent coefficients appropriate for a given range in N n (see 
Table 2). We adopted a preliminary model for this paper, but 
use of any of their final models does not affect our conclu- 
sions. 

Just as absorption due to intervening clouds will strongly 
alter the background radiation spectrum, the reemission of 
ionizing radiation from the clouds will also affect the radia- 
tion (Haardt & Madau 1996; Fardal etal 1997). Recombi- 
nations in the clouds produce continuum radiation (H I Lyc, 
He II Lyc, He II Bac, and He II 2-photon) as well as He II 


As stated previously, we neglect the possibility that the 
absorbers consist of more than one thermal and ionization 
phase to concentrate on the effects of the shape of the pho- 
toionizing radiation spectrum. We model the Lya forest ab- 
sorbers using the photoionization code CLOUDY version 90 
(Ferland 1996). The absorber clouds are modelled as plane- 
parallel slabs of constant density, illuminated on both sides 
by the ionizing radiation field. We adopt a fiducial column 
density N(H I) - 10 15 cm -2 and assume a metallicity Z 
= 10 _2 Z o . We adopt the spectral shapes for the incident 
radiation discussed in Sec. 2.7 and vary the strength of the 
incident radiation through the ionization parameter, U 
= n y ln H , where n y and n H are the number densities of ion- 
izing photons and hydrogen respectively. Each curve in Figs. 
2, 4, and 6 represents a set of photoionization models paired 
with the corresponding spectrum shown in Figs. 1, 3, and 5, 
respectively. 

Although we focus this paper on the effect of the spectral 
shape, we have also explored the effect of varying the pa- 
rameters assumed above. For the spectral shapes considered 
here, there is little change in the ratios (C ii/C IV is reduced 
by ~ 10% and Si iv/C iv is raised ~ 10%) if 7V(H I) is in- 
creased by a factor of 10, as long as the temperatures for the 
clouds are assumed to be consistent with the thermal equi- 
librium solution of the CLOUDY model. An increase (or 
decrease) by a factor of 3 in the assumed metallicity has a 
slightly larger effect. The C ii/C iv fraction is reduced 20% - 
30% with a factor of 3 reduction in metallicity and raised 
8%-15% by a factor of 3 increase in metallicity. The ratio 
Si iv/C IV is raised as much as 10% and reduced by 6% or 
less with the same changes in metallicity. Constraints on the 
H I columns of the clouds preclude large reductions in the 
assumed density for one-phase photoionization models, to 
avoid excessive length scales for the clouds and to ensure 
that contributions to the baryon density from these N(H i) 
^ 10 15 cm -2 absorbers do not exceed BBN limits. 

We normally allow the temperature to be that solved for 
by CLOUDY, assuming thermal equilibrium between photo- 
electric heating and radiative cooling. If the clouds are actu- 
ally cooler, for example, if they have expended some energy 
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Fig. 1. Mean intensity, in relative units, vs energy in Ryd. All curves have 
been normalized to the same values at v ~ ^(1 Ryd). Solid curve assumes 
IGM-filtered sources with AGN (v" 1 *) spectrum at z = 3.4. Dotted curve 
assumes a mean intensity with v' 18 above v — (1 Ryd) and v ~ v Hel , (4 

Ryd), with decrease by a factor of 2 above designed to give same 

ratio, J( 4 Ryd)/./( ! Ryd), as the properly tillered case. Short-dashed curve 
assumes an IGM-filtered spectrum: 2/3 starburst (from Sutherland & Shull 
1997) and 1/3 AGN. Long-dashed curve assumes v ' 1,8 above v = u u and 
v - v H ,. u . and a decrease by a factor of 12 above v Mell . to give same ratio, 
J{ 4 Ryd)/./< I Ryd), as the mixed-source spectrum above. 

in expansion, or if they retain thermal memory of earlier 
epochs when Compton cooling off the microwave back- 
ground was important, the ratios Si iv/C IV vs C ll/C IV may 
be increased by as much as a factor of 2 at low values ot 
C Ii/C iv (high values of U). If, as suggested by Haehnelt 
et al. (1997), the temperatures of the clouds are larger than 
that expected from photoionization thermal equilibrium, the 
Si iv fraction is more temperature-sensitive than the C iv 
fraction, and Si iv/C iv is decreased further relative to C 11/ 
C iv. These higher temperatures, suggested as a way to in- 
crease the apparent thickness of the clouds, also imply much 
smaller metallicities, as the C iv traction is decreased much 
less compared to the H I fraction by higher temperatures. For 
example, adopting the Jq ( spectrum and assuming U 
- 10" we find that increasing the cloud temperature from 
25 000 to 70 000 K decreases the C iv fraction by less than 
20% but reduces the H I fraction (and inferred metallicity) by 
a factor of 8. 

There is a further effect if the assumed temperature of the 
cloud is increased to the point that collisional ionization ot 
H 1 becomes important. In this case, He ll/H 1 may be larger 
than that assuming pure photoionization only, N( He II) may 
exceed l/a Hcll , and the cloud may become self-shielding. 
Even if this is not the case for a cloud with N( Hi) 
= 10 1:> cm” : , this becomes increasingly likely for greater 
values of N{ H I), so that radiative transfer within the indi- 
vidual clouds becomes important. In Fig. 7 we show the 
results of models that explore some of these temperature ef- 
fects but defer a complete discussion to a later paper. 

4. RESULTS 

4 A The Meta galactic Radiation Background 

Figure 2 summarizes the results of our photoionization 
models which assume different sources o! photoionizing ra- 



Fig. 2. Ratios Si iv/C iv vs C ll/C iv in photoionized models with , V { H I ) 
= K) 15 cm \ Curves match the spectra of Fig. 1. Triangles are from mea- 
sured column densities of Si iv, Civ, and C n (Songaila & Cowie 1996a. 
1996b). Squares are measurements and upper limits from Songaila & Cowie 
(1996a, 1996b). Diamonds arc measurements and upper limits from Sava- 
glio et al. (1997). Same data are also used in Figs. 4. 6, 7, 8 with model 
calculations of Si iv/C IV vs C ll/C iv. 


diation for the metagalactic background. The solid and dot- 
ted curves both assume photoionization by sources that pos- 
sess the same power-law spectral shape ( ot s = 1.8) but differ 
in the treatment of the filtering of their radiation by interven- 
ing clouds (see Fig, 1). Models incorporating a pure- AGN 
spectrum, J Q f . indicate that the Si iv/C iv ratio may be en- 
hanced by almost a factor of 2 over the simpler model for the 
absorbers with C ll/C IV > 0.1. Even in this regime. Si is still 
overabundant relative to C, as previously argued by Songaila 
& Cowie (1996a). However, a factor of two enhancement 
may be sufficient. If the absorbers are photoionized, an in- 
creasing C ll/C iv ratio is directly related to a decreasing ion- 
ization parameter (Fig. 2). At high C Ii/C iv (low U) there 
are few high-energy photons available to ionize the higher 
ionization stages of silicon. As a result, the enhanced ioniza- 
tion of Si in by He II Lya photons emitted by clouds is 
enough to increase the ratio of Si iv/C iv over that of the 
broken power-law models. Compared to more realistic cal- 
culations with cloud reemission, models for the spectral 
shape of the metagalactic background with a simple break at 
4 Ryd underestimate the number of high-energy photons in 
the metagalactic background. At low C Ii/C iv (high ( 7 ), the 
Si iv fraction is strongly depleted by the increased amount of 
high-energy photons in the realistic spectrum. 

A metagalactic ionizing background dominated by stellar 
sources, y MB = 7 S( ; / , increases the Si iv/C IV ratios at 
C n/C iv > 0.1 by 20% -50% over that of the 7 N1B = J QJ case. 
Once again, using a more realistic spectrum increases the 
Si iv/C IV ratios by almost a factor of 2 compared to spectra 
with simple attenuation of He II ionizing photons. As in the 
pure AGN case, the simple representation overestimates the 
corresponding intensity at higher energies. To be compatible 
with observed Si iv/C iv ratios if C n/C iv >0.1, the relative 
abundance of Si must be increased by approximately a factor 
of 2. In all cases, the large Si iv/C iv ratios in absorbers with 
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Fig. 3. Mean intensity, in relative units, vs energy in Ryd. All curves have 
been normalized to the same values at v - y H . All curves are mixtures 
where J MB is the IGM-filtered AGN spectrum shown Fig. 1, and J L is the 
unfiltered (“local”) stellar spectrum (Sutherland & Shull 1997). Curves 
assume ~ 0 (solid), 5 (dotted), 10 (short-dashed), and 20 (long- 

dashed). Dot-long-dashed curve assumes only local stellar radiation, J L . 

C Il/C IV < 0. 1 are difficult to explain in the context of single- 
phase models unless Si/C is enhanced by an order of magni- 
tude. Chemical evolution models with Si/C^3 (Si/C) G are 
unrealistic, even if the nucleosynthesis is dominated by mas- 
sive stars (Woosley & Weaver 1995). 

4.2 A Photoionizing Background Dominated by Local 
Sources 

One resolution to the puzzle of large Si iv/C IV at small 
C IL/C IV may be that the ionizing radiation incident on the 
absorbers is dominated by stellar sources of radiation. As 
Fig. 4 shows, if the ionizing spectrum is entirely stellar, very 
large ratios of Si iv/C iv are always possible due to the steep 
dropoff in radiation more energetic than 45 eV. In practice, 
however, as Figs. 3 and 4 show, if the metagalactic back- 



Fig. 4. Ratios of Si iv/C iv vs C N/C iv in photoionized models with 
/V(H i) = 10 |S cm Curves match the spectra of Fig. 3. 



Fig. 5. Mean intensity, in relative units, vs energy in Ryd, where all curves 
have been normalized to the same values at v=v H . All curves are mixtures 
where J MB is the filtered stellar/AGN spectrum shown in Fig. I, and J L is 
the unfiltered stellar spectrum from Sutherland & Shull (1997). Curves as- 
sume J L U MB = 0 (solid), 3 (dotted), 5 (short-dashed), 10 (long-dashed), and 
20 (dot-dashed). 

ground is primarily due to quasars (Jmb = Jqj), even if 
JlU mb = 20, high Si iv/C IV ratios at low C n/C iv require 
Si enhancements of a factor of 10. 

If the metagalactic background is dominated by stellar 
sources (Jmb = Jsqj)i large enhancements of Si may be un- 
necessary (see Figs. 5 and 6). For many absorbers with 
Cll/Civ>0.1, no overabundance of Si is necessary if 
J mb” 3 — 10. At low C Ii/c iv, J l !Jmb must exceed 10 if 
enhancements in Si exceeding 2 are to be avoided. As our 
order of magnitude estimates in Sec. 2.1 indicate, this would 
require that the absorbers lie within 30 kpc of starburst gal- 
axies. However, this situation is quite feasible, given realistic 
constraints on absorber size and metal transport distances. 

4.3 Nonequilibrium Temperatures and Nonuniform 
Radiation Fields 

As Figs. 7(a) and 7(b) show, if the temperatures of the 
clouds are cooler than the photoionization thermal equilib- 



C Il/C IV 

Fig. 6. Ratios Si iv/C iv vs C Il/C iv in photoionized models with N(H i) 
= 10 15 cm 2 . Curves match the spectra of Fig. 5. Dot-long-dashed curve 
assumes only local stellar radiation, J L . 
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C II/C IV 

Fig. 7. (a) Effects of varying cloud temperature. All cloud models assume a 
metagalaetic background with quasar sources only, * Solid curve 

adopts the equilibrium temperature computed by CLOUDY for a given ion- 
ization parameter and spectral shape (varies between T- 16 500 K and 
F-45 000 K.) Dotted curve assumes constant T= 15 (XX) K for all models; 
dashed curve assumes constant F==50 000 K for all models, (b) Same as (a), 
except background includes both quasars and hot stars, J %\n = ^ sqj' ' 

rium temperatures (Ferrara & Giallongo 1997; Zhang el al. 
1997), much higher ratios of Si lv/C IV are possible for low 
C II/C IV (high U) clouds. These high U clouds are precisely 
the clouds for which the photoionization thermal equilibrium 
temperature is highest ( T » 40 000 K). The models denoted 
by the dotted curves in Figs. 7(a) and 7(b) assume T 
= 15 000 K, which may not be compatible with the observed 
linewidths, and which are cooler than is usually assumed for 
Lya forest clouds. If higher temperatures are associated with 
the clouds, for example the dashed curves in Figs. 7(a) and 
7(b), which assume T = 50 000 K, it is more difficult to ac- 
count for the ratios in a one-phase model. 

Another way to enhance Si iv/C iv at low C II/C IV, as 
suggested by Songaila & Cowie (1996a) and Savaglio et al 
(1997), is to assume that no photons above 4 Ryd are present 
in the incident ionizing spectrum. This is the limiting case of 
a cloud embedded in a region where He n I-fronts have not 
yet overlapped. The Si iv/C iv ratios in Fig. 8 are not as high 
as those for the limiting case when all incident radiation is 
starburst radiation, since radiation between 45 and 54 eV is 
included, which will ionize Si iv to Si v. This effect is seen 
in comparing models which assume pure AGN and starburst/ 
AGN sources in Fig. 8. As Fig. 8 also shows, a careful treat- 
ment of the propagation of these He it I-fronts is necessary, 
particularly because most AGN possess x-ray emission with 
spectra that flatten to approximately v 1 above 0.3 keV. We 
have shown that including high-energy radiation with E 
> 25 Ryd is sufficient to reduce Si iv/C iv ratios substan- 
tially at low C Il/C iv. 

5. DISCUSSION 

From the results of our one-phase photoionization models 
we draw the following conclusions. 



Fig. 8. Raiios Si iv/C iv vs C n/C iv in an IGM in which He n ionization 
fronts have not overlapped. Solid curve repeats models lor the (iltcred 
starburst/AGN spectrum in Fig. 2. Dotted assumes filtered starbursi/AGN 
spectrum with all photons with E>4 Ryd removed. Long-dashed curve as- 
sumes v 1,8 AGN spectrum with all photons with E>4 Ryd removed. 
Short-dashed curve assumes filtered starburst/AGN spectrum with 
4<E<25 Ryd photons removed. 

(1) Ionizing Radiation Field . — The metagalaetic radiation 
field is likely to include both power-law (AGN) and stellar 
(hot-star) components. The Lya forest absorbers with metal 
lines may also experience a local radiation field from star- 
burst galaxies within 50-100 kpc. 

(2) Si/C Overabundance. — For plausible mixtures of stel- 
lar and AGN spectra in the metagalaetic background, our 
photoionization models produce enhanced Si iv/C iv ratios, 
consistent with high-z absorbers with Si/C fia 2(Si/0)o at low 
ionization (Cn/Civ > 0.1). For absorbers with C n/C iv 
< 0.05, it is difficult to account for the high values of Si IV/ 
Civ unless Si/C>10(Si/C) G , an unrealistically large value 
for massive-star nucleosynthesis. These systems may include 
photoionization from local stellar sources as well as hot, col- 
lisionally ionized gas. 

(3) Local Ionizing Sources. — If the radiation field incident 
on the absorbers is dominated by a nearby starburst galaxy, 
the Si iv/C IV ratios are further enhanced. If C ll/C IV > 0. 1 , 
no Si/C overabundance is necessary to explain Si iv/C IV ra- 
tios if absorbers are within about 40 kpc of a starburst galaxy 
and the background is dominated by stellar sources. Even il 
the metagalaetic background is dominated by AGNs, close 
proximity to a starburst galaxy may reduce the needed Si/C 
overabundance to a factor 1.5. When C ll/C iv<().()5, it re- 
mains very difficult to account for the high values of Si iv/ 
C IV with photoionized models, although models with locally 
dominated radiation may only require Si/C enhancements of 
a factor of 2-3. 

(4) Non-Overlapping He // I-Fronts. — Around quasars 
whose ionizing continua have steep spectral indices ( ot s 
> 1.84), the He il I-fronts will lie within the H I I-fronts. If 
many absorbers lie in regions where all photons above 4 Ryd 
are attenuated (Songaila & Cowie 1996a), we obtain good 
agreement with almost all measurements of Si iv/C IV if Si/C 
is enhanced by a factor of 2-3. However, including even a 
small contribution of higher energy photons in the back- 
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ground increases the needed Si/C overabundance to an order 
of magnitude for absorbers with the lowest C Il/C IV ratios. 

(5) Temperature Effects . — If the absorbers are cooler than 
expected for thermal equilibrium between photoelectric heat- 
ing and radiative cooling, the Si iv/C IV ratio is increased for 
a given C ll/C iv ratio. At low C il/C iv, this enhances Si IV/ 
C iv by a factor of 5 in our photoionization models. Con- 
versely, if the absorbers are hotter. Si iv/C iv is lower for a 
given C Il/C IV ratio. 

While we have described several processes that might in- 
crease the Si iv/C IV ratio, an important additional constraint 
on these possibilities is the fact that the high ratios are pref- 
erentially found at z > 3.1. Recent observations (Boksenberg 
1997) challenge this interpretation, by finding high Si iv/C IV 
ratios in absorbers at z = 2-3. If, however. Si iv/C iv rises 
above z > 3.1, this argues for a time dependence to whatever 
process increases this ratio. It is this property, as well as the 
increased Hen absorption toward Q0302-003 at z = 3.28 
(Jakobsen et ai 1994), that Songaila & Cowie (1996a) use to 
support their suggestion that the absorbers lie in regions 
where He u I-fronts have not overlapped at z > 3.1. From a 
higher resolution HSTIGHRS spectrum of Q0302-003, 
Hogan et ai (1997) find evidence for residual transmitted 
flux below the He II edge. Their 95% confidence upper limit, 
r He j, ^ 3, makes it less necessary to propose that Q0302— 003 
lies in a region where He II I-fronts have not overlapped, but 
does not preclude the suggestion that such regions existed at 
z > 3.1. If this is the case, radiation outside of He III regions 
would largely limit the range in ionization stages for carbon 
to II-IV, and for silicon to II-V. In AGN, however, the large 
fluxes of unattenuated photons above 45 eV can ionize Si iv 
and C in, so that high observed Si iv/C IV ratios would be 
difficult to explain for large V (small C Il/C iv). As a result, 
the spectrum below 4 Ryd must be dominated by stellar 
sources. This may be difficult to achieve with the known 
quasar luminosity functions at z > 3.5. In addition, soft x-ray 
radiation from quasars, which is less attenuated by interven- 
ing clouds even in the low He III porosity case, may again 
make high Si iv/C IV, low C Il/C iv absorbers difficult to un- 
derstand if they are solely photoionized by a metagalactic 
background. 

New or improved measurements of He II absorption may 
soon indicate whether He u I-fronts have overlapped well 


before z = 3. 1 . In that case, a sharp change in the shape of the 
metagalactic radiation field may be less plausible. Still, in 
general, the trend is likely to be more attenuation of photons 
above the 4 Ryd limit with increasing redshift. Other time 
dependent effects that enhance the ratio of Si iv/C iv are 
possible. Higher abundances of Si relative to C are associ- 
ated with metal yields from the most massive stars. How- 
ever, in a flat universe with 6 = 0.75, the age of the universe 
exceeds 1G 9 years by z = 3.1, time enough for lower-mass 
stars to enrich the gas with carbon. For example, if multiple 
supernovae eject metal-enriched gas into the IGM, the addi- 
tional 4X I0 8 years between z = 3.5 and z — 2.5 may make it 
more likely that the gas surrounding the high-mass stars has 
been enriched with carbon from a previous episode of star 
formation. 

The temperature effects we discuss in Sec. 4 may also 
have a time dependence. If the IGM has been uniformly 
enriched by a much earlier episode of Population III star 
formation, and if the Lya forest clouds with metals have 
their origin in growing overdensities in the IGM, then the 
absorbers observed at higher z may have formed at an earlier 
epoch. If Compton cooling off the cosmic microwave back- 
ground was a dominant coolant at this epoch (Miralda- 
Escude & Rees 1994), the cloud may have retained memory 
of this lower temperature. This effect is probably not large, 
since there is not a significant evolution in the observed lin- 
ewidths. 

Alternatively, one may resolve the puzzle of the high 
Si iv/C iv, low C Il/C IV absorbers by relaxing the single- 
phase model which we used to explore the effects of radia- 
tive transfer. As we have emphasized in a previous paper 
(Giroux et ai 1994), the production of heavy elements in 
QSO absorption systems is naturally accompanied by hot gas 
due to supernovae and hot-star winds. 

We thank Mark Fardal for the calculated metagalactic 
background radiation spectra. We also thank G. Ferland for 
discussions about updated versions of CLOUDY, and L. 
Cowie and A. Songaila for providing a summary of measure- 
ments and upper limits on metal column densities in advance 
of publication. This work was supported by the Astrophysi- 
cal Theory Program at the University of Colorado (NASA 
Grant No. NAGW-766). 
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